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Preface 


This  independent  study  began  as  an  effort  to  perform  a 
worst-case  analysis  on  the  amount  of  radioactive  fallout 
exposure  to  aircrew  members  when  flying  through  a  descending 
fallout  cloud  generated  by  a  nuclear  surface  burst.  This 
problem  has  recently  generated  a  great  amount  of  interest  by 
military  planners  during  this  time  of  debate  about  the  basing 
of  the  MX  missiles.  During  a  war,  an  attack  against  these 
sites  by  many  megaton-size  nuclear  weapons  could  present  a 
problem  to  aircrews  aboard  airborne  command  posts  and  similar 
strategic  aircraft.  Since  these  aircraft  require  extended 
periods  of  airborne  time,  penetration  of  radioactive  fallout 
clouds  generated  by  these  bursts  becomes  almost  inevitable. 

Presented  within  are  the  results  of  my  efforts  to  solve 
this  problem  by  developing  a  computer  code  for  calculating 
the  ionizing  dose  rate  of  a  radioactive  dust  cloud  as  a 
function  of  time,  and  also  the  external  dose  that  an  aircrew 
receives  during  and  after  the  cloud  transit.  A  simple  and 
efficient  procedure  for  handling  multiple  bursts  is  also 
presented. 

I  am  especially  grateful  to  Dr.  Charles  J.  Bridgman  for 
his  direction  and  inexhaustible  patience  during  this  research. 
I  am  particularly  indebted  to  my  lovely  wife,  Xathi,  and  to  my 
children  for  the  great  amount  of  patience  and  love  they  showed 
me  throughout  the  period  of  this  work.  Finally,  thanks  are 
extended  to  Ms.  Sharon  Gabriel  for  typing  this  report. 

Burl  E.  Hickman 
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Abstract 

This  report  will  evaluate  the  threat  of  radioactive 
fallout  to  which  aircrew  members  will  be  exposed  when 
flying  through  a  descending  fallout  cloud. 

A  computer  program  is  developed  for  calculating  the 
ionizing  dose  rate  of  a  radioactive  dust  cloud  as  a 
function  of  time,  and  also  the  dose  that  an  aircrew  receives 
when  flying  through  the  respective  cloud.  A  cloud  model  that 
is  patterned  after  the  AFIT  fallout  smearing  code  was  devel¬ 
oped.  A  comparison  is  made  between  the  activities  at 
various  altitudes  from  305  meters  to  9150  meters  to  provide 

r  1 

information  for  possible  re-direction  of  flight. 

The  external  ionizing  dose  to  the  aircrew  is  computed 
*  by  the  new  code  considering  the  cloud  size,  the  aircraft's 
transit  time  through  the  cloud,  and  the  ingestion  rate  of 
radioactive  particles  into  the  aircraft's  cabin.  Information 
is  also  provided  to  indicate  the  method  by  which  doses  can 
be  computed  from  a  cloud  of  multiple  bursts.  The  results 
demonstrate  that  total  dose  to  each  aircrew  member  is 
approximately  8  rems  after  flying  through  a  fallout  cloud 
one  hour  after  cloud  stabilization  of  a  1  Mt  burst,  with 
the  mission  continuing  for  eight  hours  subsequent  to  the 
cloud  transit. 

fv 

« 
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AIRCREW  IONIZING  DOSES 

FROM  RADIOACTIVE  DUST  CLOUD  GENERATED  BY  NUCLEAR  BURST 

I .  Introduction 

Background 

The  basing  site  of  the  MX  missile  system  would  be  a 
prime  target  for  an  enemy  first-strike  attack.  Aided  by 
the  prevailing  westerly  winds,  the  numerous  radioactive 
dust  clouds  generated  by  such  an  attack  would  cause  vast 
areas  of  the  United  States  to  be  exposed  to  a  considerable 
amount  of  radioactive  fallout.  Airborne  aircraft  will  be 
subjected  to  accidental  or  deliberate  penetration  of  these 
clouds  when  performing  their  wartime  missions.  Consequently, 
the  Strategic  Air  Command  (SAC)  has  expressed  a  concern 
about  the  amount  of  fallout  that  the  Airborne  Command  Post 
personnel  wculd  be  exposed  to  when  flying  through  the 
descending  fallout  particles  generated  by  these  bursts 
(Ref  16).  Real-time  information  on  the  location,  time, 
yield,  height  of  burst,  and  fallout  patterns  with  dose 
rates  would  provide  the  Airborne  Command  Post  battle  manage¬ 
ment  with  information  on  patterns  of  attack  and  on  areas 
which  would  be  hazardous  to  human  life  and  electronic 
equipment.  It  would  provide  information  necessary  for  the 
rational  direction  of  aircraft  during  flight  and  recovery, 
as  well  as  provide  insight  into  ground  facilities  which 
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may  be  unavailable  as  a  result  of  excessive  fallout.  For 
example,  a  simple  altitude  change  during  flight  might 
take  the  aircrew  from  an  extremely  hazardous  environment 
to  one  with  lesser  hazards.  Effective  aircraft  operating 
conditions  must  also  be  considered  in  such  a  decision, 
however . 

Aircraft  penetration  of  radioactive  clouds  poses 
potential  hazards  in  three  ways.  First,  the  crew  of  the 
aircraft  is  exposed  to  ionizing  radiation  through  the  air¬ 
craft's  ventilation  system  and  directly  through  the  skin 
of  the  aircraft.  Second,  the  electronic  equipment  could 
suffer  functional  damage  or  operational  upset  when  exposed 
to  ionizing  radiation  (Ref  1).  Third,  the  dust  into  the 
*  -  engine  and  aircraft  ventilation  system  could  cause  erosion 

of  engine  parts  (Ref  2),  This  study  will  focus  only  on 
the  first  hazard. 

Problem 

Flights  through  clouds  may  have  taken  place  previously, 
but  results  from  these  flights  could  not  be  found.  The 
problem  addressed  in  this  study  is  the  determination  of 
dose  rates  to  which  the  Airborne  Command  Post  personnel  would 
be  exposed  when  flying  through  the  descending  fallout  parti¬ 
cles  generated  by  nuclear  bursts.  Initially,  the  particle 
distributions  within  a  nuclear  cloud  as  a  function  of  yield, 
altitude,  and  settling  time  must  be  determined.  Onv.e  the 
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particle  distribution  is  known,  the  total  particle 
activities  can  be  used  to  compute  an  average  dose  rate 
and  the  total  dcse  to  an  aircrew  member  inside  the  aircraft’s 
cabin. 

Scope 

This  study  has  focused  on  modeling  the  potentially 
hazardous  cloud  and  considered  the  activity  that  may  enter 
and  settle  inside  of  the  aircraft.  The  fraction  of  enter¬ 
ing  particles  which  settle  out  is  a  difficult  question,  but 
it  directly  affects  the  computation  of  personnel  and  equip¬ 
ment  exposure.  Therefore,  for  simplification  purposes, 
only  the  worst-case  concept,  i.e.,  all  of  the  radioactive 
particles  that  enter  the  aircraft  as  it  passes  through  the 
cloud  remain  inside  the  cabin,  was  selected  to  illustrate 
the  upper  bound  of  the  problem.  This  concept  would  also 
indicate  whether  or  not  a  real  problem  does  exist.  Dust 
accumulation  in  the  interior  of  the  aircraft  is  primarily 
due  to  cabin  pressuriz? tion  and  air  conditioning,  and  to 
heavy-duty  electronic  cooling.  Consequently,  the  aircrews 
are  exposed  to  contamination  from  direct  skin  contact  and 
from  inhalation.  This  study  will  focus  only  on  the  whole 
body  exposure  due  to  external  gamma-ray  radiation. 

Assumptions 

Several  other  explicit  assumptions  are  made  in  this 
work.  They  are: 
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1.  The  initial  conditions  for  the  stabilized  cloud 
position  as  treated  by  Bridgman  and  Bigelow 
(Ref  3)  are  applicable. 

2.  The  activity-size  distribution  of  fallout 
particles  as  presented  by  DELFIC  (Ref  4)  is 
applicable . 

3.  No  self-shielding  from  gamma  radiation  exists 
within  the  aircraft. 

4.  There  is  no  significant  adherence  of  radio¬ 
active  dust  particles  to  the  external  skin 
of  the  aircraft. 

5.  Megaton-size  yields  will  be  used  against  the 
hardened  basing  sites. 

The  assumptions  will  be  discussed  in  greater  detail  later 
in  the  text. 

Approach 

The  mathematical  development  and  results  of  this 
analysis  of  the  radioactive  dust  cloud  model  are  presented 
in  Chapter  II.  The  results  are  presented  in  terms  of 
activity  density  in  curies  per  cubic  meter  versus  altitude 
at  a  specified  time.  The  mathematical  development  of  the 
external  dose  to  the  aircrew  is  presented  in  Chapter  III. 
The  results  are  then  presented  in  tabular  form  and  include 
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the  most  prominent  particle-size  group  that  affects 
the  aircrew  at  each  altitude  Chapter  IV  presents  the 
analysis  of  multiple  bursts,  and  Chapter  V  presents  the 
conclusions  and  recommendations  of  this  report. 


II.  Cloud  Model 


Background 

Several  different  models  have  previously  been 
developed  to  predict  the  radioactive  fallout  cloud, 
some  of  which  were  analyzed  by  Polan  in  1962  (Ref  5). 

All  of  these  models,  however,  concerned  themselves 
primarily  with  the  predicted  footprint  that  the  radio¬ 
active  fallout  makes  over  the  ground  as  the  cloud 
descends.  Of  these,  the  WSEG-10  and  the  DELFIC  models 
appear  to  have  been  the  most  popular  (Ref  6);  WSEG-10 
for  its  simplicity  and  ease  of  computation,  and  DELFIC 
for  its  predictive  capability.  Although  WSEG-10  is 
easy  to  use,  it  fails  to  account  for  fractionation,  to 
allow  for  variations  in  the  activity  particle-size 
distribution  of  the  fallout,  and  to  use  realistic 
settling  rates.  DELFIC  includes  all  of  these  provisions, 
but  the  program  is  very  expensive  to  run  and  the  model 
is  quite  difficult  to  learn. 

Patrick  (Ref  7)  and  Patrick  et  al  (Ref  2)  investi¬ 
gated  aircraft  penetration  of  radioactive  clouds  by 
using  a  simplified  cloud  model,  but  the  model  had  some 
major  shortcomings.  Specifically,  they  assumed  the 
cloud  to  have  no  variation  in  its  parameters  with  alti¬ 
tude  or  range  from  the  point  of  detonation,  with  the 
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exception  that  it  was  a  descending  volume.  Radioactive 
dust  was  depleted  from  this  fixed  cloud  by  an  effect 
called  the  "nibbling  rxiuse''  assumption,  i.e.  ,  particles 
with  the  largest  radius  fall  out  before  any  of  the 

particles  with  the  next  smaller  radius.  The  settling 

v 

rates  of  these  particles  were  computed  by  Stokes  law 
with  a  constant  air  viscosity.  Fractionation  was  also 
treated  in  a  somewhat  arbitrary  manner. 

Since  that  time,  Bridgman  and  Bigelow  have  developed 
a  new  fallout  prediction  model,  the  AFIT  model,  that 
resembles  WSEG-10  in  some  respects,  but  which  produces 
results  that  compare  better  with  those  of  DELFIC  (Ref  3). 
Unlike  WSEG-10,  this  new  model  does  account  for  variation 
in  activity  size,  fractionation,  and  uses  realistic 
settling  rates.  Like  WSEG-10,  however,  it  is  a  cloud 
smearing  model  that  concerns  itself  with  the  predicted 
fallout  footprints  over  the  ground. 

This  report  is  concerned  only  with  the  cloud  when 
it  is  airborne,  not  when  it  is  smeared  over  the  ground. 
Therefore,  the  familiar  g(t)  function,  which  is  the 
normalized  rate  of  arrival  of  activity  on  the  ground 
everywhere  and  which  is  used  in  both  the  WSEG-10  and 
AFIT  models,  was  inapplicable.  Consequently,  a  new 
model  which  is  patterned  after  the  AFIT  model  because 
of  its  advantages  over  WSEG-10  was  developed. 


Theory 

The  DELFIC  model  is  considered  to  be  the  standard- setting 
model.  Since  the  cloud  model  suggested  by  Bridgman  and 
Bigelow  compares  closely  to  the  DELFIC  model  and  computes 
in  seconds  or  less*  its  initial  conditions  were  used.  Speci¬ 
fically*  it  is  assumed  that  the  radioactive  dust  cloud  is 
generated  from  a  surface  burst.  For  this  thesis*  the  burst 
is  further  assumed  to  be  of  a  megaton-size  yield.  The  cloud 
is  formed  from  the  vaporization  of  surface  material  and  the 
weapon  itself.  The  cloud  then  rises  and  stabilizes  at  a 
maximum  height,  Hc  ,  approximately  10  minutes  after  detona¬ 
tion.  The  cloud  center  height  as  presented  by  WSEG-10  is 

H  -  44. 0+6.1 in (YLD)  -  0 . 205 (£n (YLD)+2 . 42) | *n (YLD)+2 . 42 | 

CD 

where  H  is  in  kilofeet  and  YLD  is  in  megatons 
c 

The  stabilized  cloud  is  modeled  as  a  right  circular 
cylinder,  as  shown  in  Figure  1.  The  radioactivity  in  the 
stabilized  cloud  is  assumed  to  be  normally  distributed  about 
Hc  in  the  vertical  and  normally  distributed  about  ground 
zero  in  the  horizontal  directions.  The  radioactivity  is  also 
a  function  of  the  particle  sizes  that  make  up  the  cloud  and 
a  function  of  time  t  ,  namely: 
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Figure  1.  Fallout  Cloud  Model 


z 


Particle 

Distribution 

Centers 


*4* - 

Downwind  (x) 


Figure  2.  Fallout  Cloud  at  Time  t 
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f(x,t)  - 


/7¥  o. 


X-V„t  2 


C3) 


f(y,t) 


/Tn  o. 


-(f)* 

y 


(4) 


and 


f  (z,t) 


✓7?  o. 


(5) 


where  x  is  the  downwind  distance ,  y  is  the  crosswind 
distance,  t  is  time,  z  is  the  vertical  height,  and  a 
(which  will  be  discussed  later  in  the  text)  is  a  time-varying 
standard  deviation  or  diffusion  coefficient.  Vx  ,  in  WSEG-10, 
is  assumed  to  be  the  constant  wind  which  rotates  uniformly 
with  altitude;  however,  since  the  aircraft  and  the  cloud 
are  suspended  within  the  same  medium,  Vx  is  taken  to  be 
the  true  airspeed  of  the  aircraft.  The  standard  deviations 
in  the  horizontal  directions  are  also  assumed  to  be  slowly 
varying  functions  of  time. 
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The  radioactivity  as  a  function  of  particle  sizes  is 
represented  by  100  discrete  particle- size  groups  which  were 
taken  from  the  DELF1C  model  (Ref  4).  This  size  distribution 
is  converted  to  an  activity  distribution  using  the  fraction¬ 
ation  treatment  of  Bridgman  and  Bigelow.  The  fall  mechanics 
of  the  size  groups  behave  according  to  the  equations  of 
McDonald  (Ref  9)  and  Davies  (Ref  10).  Therefore,  at  some 
later  time  t  ,  the  cloud  has  exp* ided  not  only  in  the  hori¬ 
zontal  directions,  but.  has  also  fallen  in  the  vertical  direc¬ 
tion  as  well.  This  results  in  each  particle-size  group 
distribution  being  located  at  different  altitudes,  as  shown 
in  Figure  2. 

In  this  report,  two  (2)  times  the  standard  deviation, 

a  ,  in  the  horizontal  and  vertical  directions  is  considered 

to  be  the  lateral  limit  from  the  center  of  the  cloud.  This 

limit,  when  taken  in  either  direction,  contains  from  95  to  98 

percent  of  the  total  activity  of  the  cloud  in  any  direction. 

Values  for  o„  and  om  have  been  derived  for  the  WSEG-10 
y  z 

and  AFIT  models  (Refs  3,  11  and  12);  however,  because  the 
g(t)  function  was  used,  a  a%  function  was  never  calculated. 
Therefore,  it  is  assumed  that  the  x  direction  will  behave 
in  a  similar  fashion  as  the  torroidal  growth  of  the  y 
direction;  but,  a  shear  term  will  not  be  required  when  assuming 
a  single  effective  wind  in  the  x  direction.  The  equation 
for  ox  then  is 
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t 


o02  Cl  ♦  8T/TC] 


(6) 


a 


where  T  is  the  time  of  arrival,  not  to  exceed  three  hours, 

and  T  is  the  time  constant, 
c 

Now  that  the  basic  equation  for  the  cloud  model  is  known, 
values  for  specific  altitudes  and  specific  times  will  be 
required  to  solve  the  original  problem.  Therefore,  the  lat¬ 
eral  limits  of  the  cloud  and  the  particle  sizes  that  make  up 
the  cloud  at  a  specified  altitude  are  required  at  the  time 
(TA)  the  aircraft  arrives  at  the  cloud.  The  f(z,t)  function, 
then,  is  also  a  function  of  particle  size.  Equation  (5)  then 

becomes  t 

za/c-<zo 
-IjC - 

f(z,t,r)  -  - i -  e 

✓2?  az 

where  z  >  is  the  aircraft  altitude,  z^  is  the  original 
a/c  *  o  ° 

cloud  center  height,  and  V  (r)  is  the  fall  velocity  as  a 

function  of  particle  size.  This  expression,  written  for  each 

of  the  100  activity-size  groups,  will  provide  the  amount  of 

activity  from  each  group  at  a  given  altitude.  Since  the  cloud 

is  of  finite  extent,  the  total  cloud  height  is  considered  to 

be  four  times  the  standard  deviation  in  the  vertical  direction. 

Consequently,  the  extent  of  the  contribution  of  each  particle- 

size  group  depends  upon  whether  2o  of  a  given  particle-size 

£ 
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group  distribution  lies  within  the  aircraft's  altitude. 

The  cloud  activity  at  its  horizontal  center  at  a  given  TA 
is 


100 


A.  t 


l  Ai 


-h 


i"l 


(8) 


where  t  is  in  hours  after  detonation,  and 


A* 


Ai  *  XW  f  (x,t)f  (y,t)f  (z,t,r) 


(9) 


A*  ,  which  is  the  activity  in  disintegrations  per  second  at 
one  hour,  is  given  by 


A*  -  530(3.7  x  10l8)YLDKT  x  FF 


(10) 


FF  is  the  fission  fraction  and  YLD  is  the  yield  in  kilo- 
tons.  Caution  must  be  exercised  at  this  point  because  all 
activities  are  unit  time  reference  activities.  They  must  be 
transferred  to  real  time  using  the  Way-Wigner  decay  formula. 

Note  that  Eq  (10)  is  the  activity  per  volume  at  the 
horizontal  center  of  the  cloud.  Figure  3  shows  the  activities 
of  a  1  Mt  burst~at““altT^de^  from  9150  meters ,  approxi¬ 

mately  30,000  feet  MSL,  to  305  meters,  approx imately~l7<KH*~ — ... 
feet  MSL,  and  at  aircraft  arrival  times  at  the  horizontal 
center  of  the  cloud  that  range  from  0.5  hours  to  eight  hours 
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DENSITY.  (CURIE  PER  CUBIC  METER) 


after  cloud  stabilization.  Notice  that  the  activity  per 
cubic  meter  decreases  at  lower  altitudes.  This  suggests 
that  a  lessor  dose  rate  would  be  received  by  the  aircrew 
if  lower  altitudes  were  selected.  These  figures*  however, 
do  not  include  any  stem  fallout.  The  values  do  show  the 
activity  that  will  be  ingested  inside  the  aircraft  at  the 
given  times  from  a  1  Mt  burst.  Considerations  for  multiple 
bursts  are  developed  in  Section  IV.  Cabin  dose  rates  will 
be  discussed  in  Section  III. 


IS 


III.  Cabin  Analysis 


Since  SAC  is  interested  in  the  exposure  to  Airborne 
Command  Post  personnel ,  an  airborne 'Command  type  aircraft, 
such  as  the  KC-135,  is  used  in  this  study.  This  aircraft 
cruises  at  231.5  m/s.  It  has  an  internal  volume  of 
246.357  m1  (Ref  14)  and  a  mass  flow  rate  of  outside  air 
into  its  air  conditioning/pressurizati^n  system  of  approxi¬ 
mately  150  lbs/min  (Ref  15).  Therefore,  as  the  aircraft 
sweeps  through  the  cloud  in  the  x  direction,  it  will  be 
ingesting  at  a  rate  of  150  lbs/min  those  particle  sizes 
that  are  co-located  at  the  aircraft’s  altitude.  The  crew 
will  be  exposed  to  this  ingested  radioactivity  and  to 
whatever  sky  shine  that  is  present. 

Sky -Shine 

As  the  aircraft  passes  through  the  radioactive  cloud, 
the  aircrew  will  be  exposed  to  direct  gamma  radiation  that 
is  passing  directly  through  the  aircraft's  skin  and  through 
any  equipment  or  material  inside  the  aircraft,  because  of 
the  zero  self -shielding  assumption.  Initially,  the  radia¬ 
tion  would  start  at  a  low  level  as  the  aircraft  enters  the 
cloud,  increase  to  a  maximum  level  as  the  aircraft  reaches 
the  horizontal  center  of  the  cloud,  and  then  decrease  again 
as  it  exits  the  other  side  of  the  cloud.  A  worst-case 
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analysis  can  be  made  if  the  analysis  was  performed  at  the 
cloud’s  horizontal  center  for  the  time  duration  that  is 
required  to  traverse  the  cloud.  Although  it  is  not  precisely 
true,  it  is  assumed  that  the  activity  of  the  dust  particles 
is  uniform  vitbin  a  few  mean  free  paths  of  the  gamma  radiation 
at  a  gi/er.  alcitule.  Therefore,  ignoring  th?  shielding  by 
the  structure  and  internal  equipment,  tne  dose  rate  can  be 
calculated  by  the  use  of  spherical  integration,  namely 


-y  s 

y  a  e  s  2sin(frded(frds 
P  4  TVS  2 


i  l) 


Pfl 

where  —  is  the  absorption  coefficient  of  tissrv?,  and  \xf 
is  the  attenuation  coefficient  of  air  and  s  is  the  radius 
of  the  cloud.  This  equation  readily  equates  tu 


1.6  x  10 


1 1 


100 


i-1 


a 


l  A.  —  (1  - 

*•_  i  mi  ' 


-V5, 


PVi 


(12) 


By  letting  s  approach  »  ,  which  again  is  a  worst-case 
assumption,  Eq  (12)  reduces  to 


1.6  x  10 


i  i 


100 
T  A. 
i-1  1  PWt 


a 


(1~ 


The  total  dose  then  can  be  computed  by  integrating  Eq  (13) 
vith  respect  to  '.ime  with  the  limits  of  TA  to  TA  *  A«  , 


17 


where  At  is  the  transit  time  through  the  cloud: 
t0+At 

D  -  /  Dt’1*2  dt'  (14) 

*a 

Tables  I  through  IX  provide  the  results  of  sky-shine  from  a 
1  Mt  burst  at  different  arrival  times.  It  can  be  seen  that 
the  total  dose  is  significant  at  early  times,  even  though 
the  length  of  time  required  to  traverse  the  cloud  is  brief. 

Ground- Shine 

t  has  been  shown  that  sky -shine  is  a  significant  factor 
iwh  reference  to  this  problem.  On  the  other  hand,  it  can 
be  shown  that  the  ground-shine  (radioactivity  from  the  cloud 
particles  that  have  fallen  to  the  ground)  is  an  insignificant 
factor,  provided  the  aircraft  is  flying  at  least  a  few  gamma 
ray  mean  free  paths  above  the  ground.  Patrick  researched 
the  problem  in  1976  (Ref  8);  however,  he  admitted  that  for 
short  flyovers,  as  is  the  case  in  this  study,  his  results 
were  questionable.  Since  his  report,  no  additional  studies 
of  this  problem  have  been  found. 

Consider  now  that  the  aircraft  is  at  height  L  above 
the  surface  radiation  as  in  Figure  4.  An  assumption  that 
consider?  the  surface  activity  to  be  uniformly  distributed 
is  made  in  order  to  simplify  the  problem.  In  reality,  the 
surface  activity  involves  a  more  complex  function  because 
the  activity  varies  with  location  and  with  time.  However, 


1C 


Side  View 


Top  View 


\ 


making  this  assumption  reduces  the  problem  to 

dCD]  -  (.M*,yWArea).  V  )  2„dr  C15) 

4tts2 

where  s2  *  r2  +  L2  and  L  is  a  constant  value.  A(x,y)  is 
the  activity  at  the  one-hour  point  of  those  particles  that 
have  reached  the  surface.  Accounting  for  units,  Eq  (15) 
reduces  further  to 

00 

D  *  A(x,y)  9.195  x  10“**  /  e^  dz  (16) 

z  z 

where  z  *  y^L  .  When  the  altitude  ranges  of  9150  meters 
to  305  meters  for  this  study  are  considered,  values  on  the 
order  of  A(x,y)  x  10"  2 9  and  A(x,y)  x  10"  1 1  are  obtained. 
Consequently,  ground-shine  is  omitted  in  the  final  analysis 
of  this  study. 

Cabin  Ingestion 

The  cloud's  activity  at  a  given  altitude  was  calculated 
in  terms  of  the  activity  per  cubic  meter.  When  the  aircraft 
traverses  the  cloud,  it  travels  completely  through  the  cloud 
in  the  x  direction,  and  this  happens  regardless  of  the 
time  after  detonation.  Therefore,  the  F(x,t)  function  in 
Eq  (2)  is  always  equal  to  unity  when  integrated  over  the 
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transit  time  TA  -  2 ax/Vx  to  TA  +  2cx/Vx  .  Then  Eq  (2) 
reduces  to  the  activity  (DPS)  per  square  meter,  namely 

A/m2  (outside  air)  *  530  x  3.7  x  101* 

x  S00F(y,t)F{z,t,r)  (17) 

m2 

The  cabin  air  Ls  taken  through  the  fifth  stage  of  the  engine 
compressor  in  such  a  manner  that  a  constant  mass  flow  of  150 
pounds  of  air  per  minute  enters  the  cabin  (Ref  15).  This 

% 

mass  rate  is  the  same  at  any  aircraft  velocity  or  altitude. 
Thus,  the  equivalent  area  of  the  effective  cabin  air  inlet 
port  is 

Inlet  Area  -  x  - 1 -  x  - ± - 

min  Pa:  */£t3  V  Ft/min 

all  a 

and  the  total  activity  ingested  per  pass  is  the  area  times 
Eq  (17);  namely 

Activity  inside  cabin  (DPS)  * 

x 

where  p  is  the  air  density  at  the  aircraft  altitude  (m) 
and  Vx  is  the  aircraft  velocity  in  meters  per  hour.  This 


A/m2 (outside  air) 


150 

“ — V 
Kair  x 


(18) 


equation  is  based  upon  a  gamma  ray  energy  of  1  Mev  (Ref  l7)* 
Note  that,  as  the  velocity  V  increases,  the  activity 

A. 

that  enters  the  cabin  decreases. 

Modeling  the  behavior  of  the  dust  particles  after 
they  have  entered  the  aircraft  presents  another  complication 
because  of  the  placement  of  air-conditioning  duct  outlets  and 
electronic/mechanical  equipment.  By  considering  the  inges¬ 
tion  flow  rate  and  the  total  internal  volume  of  the  cabin, 
it  was  determined  that  the  absolute  minimum  time  for  S9I  of 
the  internal  air  to  be  displaced  is  approximately  21  minutes. 
Therefore,  it  is  assumed  that  the  activity  that  enters  the 
cabin  is  equally  distributed.  The  cabin  activity  per  cubic 
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meter  is  then  the  total  activity  divided  by  the  total 
volume  of  246.357  cubic  meters. 

To  obtain  meaningful  dose  rates,  the  internal  cabin  is 
modeled  as  a  right  circular  cylinder  having  the  dimensions 
of  the  internal  cabin.  Again,  self -shielding  is  considered 
to  be  zero,  and  the  assumption  that  all  the  activity  that 
enters  the  cockpit  remains  inside  the  cabin  is  still  valid. 
The  dose  rate  at  the  center  of  the  cylinder  due  to  the  S’*  - 
pended  dust  is 


D 


16.0274 
2  / 

0 


2.7432  2  it 

/  /  A/ 

0  6  pa 


-Ut  (r^+z2)1* 

e - rdOdrdz  (It; 

4ir  (r2+z2) 


where  r  is  the  radius  of  the  cylinder,  z  is  one -half  the 
height  of  the  cylinder,  and  A  is  the  one-hour  activity  of 
those  particles  that  enter  the  aircraft.  Due  to  the  complex 
nature  of  the  integral  in  Eq  (19),  a  numerical  integration 
is  required.  After  integrating  over  de  ,  the  equation 
becomes 


16.0274 


/ 

0 


2.7432  -Vt(r2+z2)* 

/  2 - 

0  (r2+z2) 


rdrdz 


(20) 
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and 


D  -  3.47  A  (21) 

tissue 

when  the  numerical  integration  is  performed.  The  total 
dose  is  obtained  from  Eq  (14).  The  dose  results  from  a 
1  Mt  burst  at  aircraft  arrival  times  ranging  from  0.5  hours 
to  8  hours  are  shown  in  Tables  I  through  IX. 
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TABLE  I 


Sunmary  of  Doses  and  Prominent  Particle  Sites 
at  1  Mt  Burst,  TA  ■  O.S  Hrs ,  8-Hour  Mission  Flown  After  Cloud  Transit 


Altitude 

(meters) 

Internal  Cabin  Dose 

(rems) 

Sky- Shine  Dose 
(rons) 

Most  Prominent  Particle 

Size  Diameter  (vm) 

9150 

4.544 

18.336 

196.38 

8840 

3.993 

16.116 

216.52 

8535 

3.499 

14.120 

228.62 

8230 

2.760 

11.135 

241.79 

7925 

2.305 

9.300 

253.68 

7620 

2.017 

8.141 

268.29 

7315 

1.770 

7.142 

283.74 

7010 

1.586 

6.399 

302.50 

6710 

1.426 

5.757 

317.37 

6400 

1.290 

5.207 

338.35 

6100 

1.182 

4.768 

338.35 

5800 

1.064 

4.295 

360.71 

5490 

0.975 

3.934 

390.75 

5180 

0.894 

3.609 

390.75 

4880 

0.815 

3.289 

416.60 

4570 

0.749 

3.022 

416.60 

4270 

0.690 

2.784 

449.83 

3960 

0.636 

2.565 

478.79 

3660 

0.587 

2.367 

478.79 

3355 

0.540 

2.179 

525.35 

3050 

0.505 

2.040 

525.35 

2745 

0.467 

1.883 

525.35 

2440 

0.431 

1.738 

573.68 

2135 

0.397 

1.600 

573.68 

1830 

0.372 

1.501 

631.48 

1525 

0.343 

1.386 

631.48 

1220 

0.323 

1.305 

631.48 

915 

0.298 

1.204 

700.69 

610 

0.280 

1.130 

700.69 

305 

0.257 

1.040 

37349.0 

aciwr  "  <  '  ■ j<  ' \ * it i  i  r<T 
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TABLE  II 


Sunmary  of  Doses  and  Prominent  Particle  Sizes 
at  1  Mt  Burst,  TA  *  1.0  Hrs,  8-Hour  Mission  Flown  After  Cloud  Transit 


Altitude 

(meters) 

Internal  Cabin  Dose 

(rents) 

Sky -Shine  Dose 
(rents) 

Most  Prominent  Particle 

Size  Diameter  (tan) 

9150 

2.360 

5.791 

121.0 

8840 

2.101 

5.156 

125.  s 

8535 

1.865 

4.578 

137.0 

8230 

1.535 

3.767 

142.6 

7925 

1.309 

3.212 

149.62 

7620 

1.160 

2.846 

157,00 

7315 

1.043 

2.559 

163.4 

7010 

0.946 

2.322 

170.05 

6710 

0.867 

2.128 

178.41 

6400 

0.788 

1.934 

187.18 

6100 

0.730 

1.792 

196.38 

5800 

0.672 

1.648 

196.38 

5490 

0.622 

1.527 

206.04 

5180 

0.577 

1.416 

216.52 

4880 

0.537 

1.317 

228.62 

4570 

0.498 

1.223 

228.62 

4270 

0.464 

1.138 

241.79 

3960 

0.431 

1.057 

253.68 

3660 

0.401 

0.985 

253.68 

3355 

0.378 

0.928 

268.29 

3050 

0.353 

0.865 

268.29 

2745 

0.328 

0.806 

283.74 

2440 

0.310 

0.761 

283.74 

2135 

0.290 

0.711 

302.50 

1830 

0.270 

0.663 

317.37 

1525 

0.256 

0.629 

317.37 

1220 

0.239 

0.588 

338.35 

915 

0.227 

0.557 

338.35 

610 

0.212 

0.520 

360.71 

305 

0.197 

0.483 

360.71 
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TABLE  III 


"1 


f.i 

jr  Suninary  of  Doses  and  Prominent  Particle  Sizes 

|  at  1  Mt  Burst,  TA  ■  2.0  Hrs,  8-Hour  Mission  Flown  After  Cloud  Transit 


L 

h 

r 

Altitude 

(meters) 

Internal  Cabin  Dose 

(rems) 

Sky-Shine  Dose 
(rems) 

Most  Prominent  Particle 
Size  Diameter  (un) 

1 

9150 

0.962 

1.525 

80.163 

i 

8840 

0.855 

1.361 

83.435 

8535 

0.771 

1.222 

86.840 

: 

8230 

0.648 

1.028 

90.383 

► 

7925 

0.SS8 

0.885 

94.07 

■ 

7620 

0.497 

0.788 

98.00 

\ 

7315 

0.451 

0.715 

101.91 

7010 

0.410 

0.651 

106.10 

6710 

0.378 

0.599 

110.40 

■ 

6400 

0.348 

0.551 

110.40 

6100 

0.324 

0.514 

116.00 

5800 

0.300 

0.476 

121.00 

1 

5490 

0.283 

0.449 

125.50 

I 

5180 

0.262 

0.416 

125.50 

4880 

0.249 

0.394 

131.43 

k 

4570 

0.231 

0.366 

137.00 

-  | 

4270 

0.219 

0.348 

142.60 

k! 

3960 

0.206 

0.326 

142.60 

b 

a 

* 

3660 

0.196 

0.311 

149.62 

t  • 1 

i  • 

3355 

0.184 

0.292 

149.62 

t  \ 

in 

3050 

0.175 

0.278 

157.00 

*  ■ 

t  ^ 

2745 

0.165 

0.261 

163.40 

b 

2440 

0.155 

0.246 

163.40 

s 

2135 

0.147 

0.234 

170.05 

1830 

0.138 

0.219 

178.41 

1525 

0.133 

0.212 

178.41 

1220 

0.125 

0.199 

187.18 

915 

0.118 

0.186 

187.18 

610 

0.112 

0.178 

195.38 

|  305 

i 

0.107 

0.169 

196.38 

I 
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TABLE  IV 


Summary  of  Doses  and  Prominent  Particle  Sizes 
at  1  Mt  Burst,  TA  ■  3.0  Hrs,  8-Hour  Mission  Flown  After  Cloud  Transit 


Altitude 

(meters) 

Internal  Cabin  Dose 

(rems) 

Sky -Shine  Dose 
(rems) 

Most  Prominent  Particle 
Size  Diameter  (um) 

9150 

0.511 

0.650 

63.06 

8840 

0.461 

0.586 

65.632 

8535 

0.413 

0.525 

68.311 

8230 

0.353 

0. 449 

71.100 

7925 

0.303 

0.385 

74.00 

7620 

0.271 

0.345 

77.02 

7315 

0.246 

0.312 

80.163 

7010 

0.223 

0.284 

83.435 

6710 

0.207 

0.263 

83.435 

6400 

0.190 

0.241 

86.84 

6100 

0.177 

0.224 

90.384 

5800 

0.165 

0.209 

94.07 

5490 

0.153 

0.195 

94.07 

5180 

0.144 

0.184 

98.00 

4880 

0.134 

0.170 

101.91 

4570 

0.125 

0.159 

101.91 

4270 

0.119 

0.151 

106.10 

3960 

0.112 

0.142 

110.40 

3660 

0.106 

0.135 

110.40 

3355 

0.100 

0.127 

116.00 

3050 

0.096 

0.122 

116.00 

2745 

0.090 

0.115 

121.00 

2440 

0.086 

0.110 

121.00 

2135 

0.082 

0.104 

125.50 

1830 

0.079 

0.101 

131.43 

1525 

0.075 

0.095 

131.43 

1220 

0.071 

0.090 

137.00 

915 

0.068 

0.087 

137.00 

610 

0.065 

0.083 

142.60 

305 

0.062 

0.078 

142.60 

TABLE  V 


Summary  of  Doses  and  Prominent  Particle  Sizes 
at  1  Mt  Burst,  TA  ■  4.0  Hrs ,  8 -Hour  Mission  Flown  After  Cloud  Transit 


Altitude 

(meters) 

Internal  Cabin  Dose 

(rems) 

Sky-Shine  Dose 
(rems) 

Most  Prominent  Particle 

Size  Diameter  (vm) 

9150 

0.324 

0.358 

53.74 

8840 

0.291 

0.322 

55.928 

8535 

0.263 

0.291 

58.211 

8230 

0.225 

0.249 

60.59 

7925 

0.197 

0.217 

63.06 

7620 

0.176 

0.195 

65.632 

7315 

0.160 

0.177 

68.311 

7010 

0.146 

0.162 

68.311 

6710 

0.133 

0.148 

71.10 

6400 

0.124 

0.137 

74.00 

6100 

0.115 

0.127 

77.02 

5800 

0.107 

0.119 

77.02 

5490 

0.100 

0.110 

80.163 

5180 

0.092 

0.102 

83.435 

4880 

0.087 

0.096 

83.435 

4570 

0.081 

0.090 

86.84 

4270 

0.077 

0.085 

90.384 

3960 

0.073 

0.081 

90.384 

3660 

0.069 

0.076 

94.07 

3355 

0.064 

0.071 

94.07 

3050 

0.061 

0.068 

98.00 

2745 

0.058 

0.064 

98.00 

2440 

0.056 

0.062 

101.91 

2135 

0.053 

0.058 

101.91 

1830 

0.050 

0.056 

106.10 

1525 

0.048 

0.053 

106.10 

1220 

0.045 

0.050 

110.40 

915 

0.044 

0.049 

110.40 

610 

0.042 

0.046 

116.00 

305 

0.039 

0.043 

116.00 
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TABLE  VI 


Summary  of  Doses  and  Prominent  Particle  Sizes 
at  1  Mt  Burst,  TA  ■  5.0  Hrs,  8-Hour  Mission  Flown  After  Cloud  Transit 


Altitude 

(meters) 

Internal  Cabin  Dose 
(rems)  x  10"2 

Sky- Shine  Dose 
(rems)  x  10"2 

Most  Prominent  Particle 

Size  Diameter  (ym) 

9150 

21.81 

21.88 

47.279 

8840 

19.75 

19.81 

49.604 

8535 

17.89 

17.95 

51.63 

8230 

15.47 

15.52 

53.74 

7925 

13.50 

13.55 

55.928 

7620 

12.01 

12.05 

58.211 

7315 

10.97 

11.00 

58.211 

7010 

9.97 

10.  OC 

60.59 

6710 

9.24 

9.28 

63.06 

6400 

8.50 

8.53 

65.632 

6100 

7.92 

7.94 

65.632 

5800 

7.39 

7.42 

68.311 

5490 

6.89 

6.91 

71.10 

5180 

6.43 

6.45 

71.10 

4880 

6.02 

6.04 

74.00 

4570 

5.69 

5.71 

77.02 

4270 

5.34 

5.35 

77.02 

3960 

5.06 

5.07 

80.163 

3660 

4.75 

4.76 

80.163 

3355 

4.51 

4.53 

83.435 

3050 

4.24 

4.25 

83.435 

2745 

3.97 

3.99 

86.84 

2440 

3.84 

3.85 

86,84 

2135 

3.61 

3.62 

90.384 

1830 

3.45 

3.46 

90.384 

1525 

3.24 

3.26 

94.07 

1220 

3.15 

3.16 

94.07 

915 

2.98 

2.99 

98.00 

610 

2.81 

2.82 

98.00 

305 

2.69 

2.70 

101.91 
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TABLE  VII 


Suninary  of  Doses  and  Prominent  Particle  Sizes 
at  1  Mt  Burst,  TA  ■  6.0  Hrs,  8 -Hour  Mission  Flown  After  Cloud  Transit 


Altitude 

(meters) 

Internal  Cabin  Dose 
(rems)  x  10"2 

Sky-Shine  Dose 
(rems)  x  10"2 

Most  Prominent  Particle 

Size  Diameter  (von) 

91  SO 

15.563 

14.522 

43.644 

8840 

14.030 

13.096 

45.430 

8535 

12.736 

11.884 

47.279 

8230 

10.995 

10. 260 

47.279 

7925 

9.646 

9.001 

49.604 

7620 

8.628 

8.051 

51.630 

7315 

7.879 

7.352 

53.740 

7010 

7.152 

6.674 

55.928 

6710 

6.636 

6.193 

55.928 

6400 

6.097 

5.689 

58.211 

6100 

5.739 

5.356 

60.590 

5800 

5.302 

4.947 

60.590 

5490 

5.002 

4.667 

63.060 

5180 

4.620 

4.311 

65.632 

4880 

4.378 

4.085 

65.632 

4570 

4.099 

3.825 

68.311 

4270 

3.889 

3.629 

68.311 

3960 

3.647 

3.403 

71.100 

3660 

3.423 

3.194 

74.000 

3355 

3.298 

3.077 

74.000 

3050 

3.097 

2.890 

77.020 

2745 

2.948 

2.751 

77.020 

2440 

2.772 

2.587 

80.163 

2135 

2.602 

2.428 

80.163 

1830 

2.523 

2.355 

83.435 

1525 

2.375 

2.217 

83.435 

1220 

2.270 

2.118 

86.840 

915 

2.175 

2.029 

86.840 

610 

2.050 

1.913 

90.384 

305 

1.927 

1.800 

90.384 

TABLE  VIII 


Sunmary  of  Doses  and  Prominent  Particle  Sizes 


at  1  Mt  Burst.  TA  »  7.0  Hrs,  8-Hour  Mission  Flown  After  Cloud  Transit 


TABLE  IX 


Summary  of  Doses  and  Piominent  Particle  Sizes 
at  1  Mt  Burst.  TA  =  8.0  Hrs ,  8-Hour  Mission  Flown  After  Cloud  Transit 


Altitude  Internal  Cabin  Dose  Sky-Shine  Dose  Most  Prominent  Particle 
(meters)  (rems)  x  10  (rems)  x  10  Size  Diameter  (pm) 


IV.  Multiple  Bursts 


The  previous  sections  of  this  study  have  been  concerned 
with  a  single  burst  only.  However,  airborne  aircraft  may  be 
subjected  to  a  radioactive  dust  cloud  generated  by  multiple 
nuclear  bursts.  The  single  burst  model  must  then  be  modified 
to  account  for  this  possibility. 

Consider  a  proposed  MX  base  hit  by  300  one-megaton 
surface  bursts.  The  burst  field  is  then  approximated  by  a 
100  Km  square  geometry  with  17.32  weapons  in  the  North-South 
(NS)  direction,  and  17.32  weapons  in  the  East-West  (EW) 
direction.  The  cloud  activity  density  will  then  be  repre¬ 
sented  as 


AVoi(x,y,z,t)  =  530  x  3.7  x  1016 


•  YLD  •  FF  •  f (x,t)f Cy,t)f (z,t)  DPS/nr 


(22) 


where,  as  Crandley  states  (Ref  13), 


f(y.t)  -  ^  I 


+  50 


100  -50  ✓Z?  a  (t) 


e  ^  dy " 


(23) 


which  is  the  sum  of  the  overlapped  gaussian  distributions. 
The  f(x,t)  function  in  the  EW  direction  is 
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Therefore,  the  description  of  the  bursts  in  the  x,y  directions 
is  a  superposition  of  300  gauss ian  distributions  laid  out  in 
a  square  matrix.  Since  the  aircraft  is  again  traveling  in 
the  x  direction,  Eq  (24)  must  be  integrated  over  the  distance, 
namely 


J  f (x,t)dx  = 


j50  i 

»  ^  -50  /?¥  ox(t) 


dx'dx 


(2  5) 


Unlike  the  single  burst  case  which  integrated  to  unity,  this 
equation  integrates  to  the  square  root  of  the  number  of  bursts, 
/300  .  Because  the  height  of  the  cloud  depends  only  upon 
yield  of  each  burst,  the  f(z,t)  function  remains  the  same 
as  a  single  burst 

100 

f (z  *t )  -  l  f.(z,t)  (26) 

i“l 

The  activity  per  square  meter  resulting  from  an  EW  pass  at 
the  horizontal  center  is 

A (DPS/m2)  «  530  x  3.7  x  101* 

•  YLD  •  FF  •  (300) (1) (l/100)f (z,t)  (27) 
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where  F(y,t)  *  at  y  ■  0  .  This  equation  is  the 

same  equation  for  a  single  burst  (Eq  (17))  of  Section  III, 
except  for  the  number  of  bursts  divided  by  field  length 
multiplied  by  the  quantity  /2?  oy(t)  .  It  can  be  concluded, 
then,  that  multiple  bursts  can  be  treated  by  multiplying  the 
results  for  a  single  burst  by  the  number  of  bursts  times 
the  standard  deviation  in  the  y  direction  and  /Fir  divided 
by  the  field  length.  In  this  example,  the  value  of  this 
constant  is 


300  bursts 
"100  Km 


/2i"  cry(t) 


where  cry(t)  is  in  kilometers. 


(28) 
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V.  Conclusions  and  Recommendations 


Conclusions 

The  results  obtained  from  this  analysis,  which  was 
done  from  a  worst-case  approach,  demonstrate  that  aircrew 
exposure  to  radioactive  fallout  generated  by  nuclear  surface 
bursts  can  be  significant.  This  significance  is  dependent 
upon  the  aircraft  arrival  time,  speed,  and  altitude.  The 
problem  is  most  prominent  at  early  arrival  and  penetration 
times  where  the  higher  sky-shine  and  the  ingestion  of  the 
more  concentrated  dust  particles  occurs.  The  hazards  can 
be  further  exaggerated  if  multiple  bursts  and/or  longer 
missions  are  flown  after  cloud  penetration.  Multiple  burst 
activities  are  determined  by  multiplying  the  results  of  a 
single  burst  by  a  constant  function. 

The  results  also  showed  that  the  radioactive  dust 
cloud’s  activity  density  decreased  with  a  decrease  in  alti¬ 
tude  at  a  given  time  after  the  burst.  When  cloud  penetration 
at  early  times  is  required,  the  heavier  particles  at  the 
lower  altitudes  are  more  widely ^distributed  than  the  lighter 
particles  at  the  higher  altitudes.  This  results  in  a  lower 
cloud  activity  density  at  the  lower  altitudes  and,  therefore, 
a  lesser  number  of  particle  sizes  would  enter  the  aircraft. 
Secondly,  filtration  of  the  larger  particle-size  groups  at 
the  lower  altitudes  is  easier  to  accomplish  than  filtration 
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of  the  smaller  particle-size  groups.  However,  filtration 
would  create  a  point  source  of  radiation  until  the  filters 
are  discarded. 

Recommendat ions 

There  are  three  recommendations  to  be  made.  First,  it 
is  recommended  that  a  stem  fallout  analysis  be  made  because 
this  analysis  was  performed  without  consideration  to  stem 
fallout.  Although  this  fallout  is  assumed  to  be  caused  pri¬ 
marily  by  larger  particle-size  groups,  the  activity  densities 
at  lower  altitudes  could  be  appreciably  affected  by  it  at 
earlier  times.  Secondly,  this  analysis  considered  a  single 
effective  westerly  wind  at  all  altitudes.  As  more  informa¬ 
tion  is  received  on  the  treatment  of  horizontal  winds,  a 
study  could  be  performed  to  consider  the  effects  on  the  cloud 
by  these  winds  at  the  changing  altitudes.  Third,  a  more 
detailed  analysis  of  dust  particle  behavior  could  be  made  to 
determine  what  percentage  of  the  radioactive  dust  particles 
actually  settled  out  inside  of  the  aircraft.  Although  a 
worst-case  analysis  was  made  in  this  study,  surely  some 
percentage  of  the  radioactive  dust  exits  the  aircraft  with 
the  outgoing  air. 
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APPENDIX  A 


Glossary  of  FORTRAN  Terms 


Input  s : 


FFR 

* 

Fission  fraction 

YLD 

- 

Yield  in  kilotons 

RT 

m 

Range  of  aircraft  from  fallout 

cloud  in  meters 

VAC 

- 

True  airspeed  of  aircraft  in  meters  per  second 

TR 

a 

Time  remaining  in  mission  after 

cloud  arrival 

Terms : 

ACALT 

s 

Aircraft  altitude 

(meters) 

HC 

ac 

Cloud  center  height 

(kilofeet) 

HC1 

- 

HC  x  3.048  x  102 

(meters) 

TA 

3 

Time  of  aircraft  arrival 

to  fallout  cloud 

(seconds) 

HR 

a 

TA/3600 

(hours) 

SZ 

« 

a2  -  .18  HC1 

(meters) 

RHOPAR 

- 

Particle  density 

(lcg/ms) 

DIA 

e 

Particle  diameters 

(meters) 

R 

s 

Radius  of  particle 

(meters) 

TC 

ac 

Time  constant 

(hours) 

SO 

m 

°o 

(meters) 

ZDIST 

'Sm 

Particle-size  distribution 

center  height 

(meters) 
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FZ 

m 

F(z,t,r) 

(meters" 1 ) 

SY 

m 

°y 

(meters) 

AR 

- 

Particle-size  group  activity 

(DPS/m3) 

AMAX 

m 

Total  activity  of  cloud 

(DPS) 

ATOT 

* 

Activity  inside  aircraft 

(DPS) 

APERM3 

- 

ATOT  per  cubic  meter 

(DPS/ms) 

ATOTTA 

- 

Cloud  activity  at  time  t 

after  detonation 

(DPS/m3) 

RHOALT 

s 

Air  density  at  altitude 

(kg/m3) 

ACPALT 

- 

Aircraft  pressure  altitude 

(Newtons /m 2 ) 

MUARHO 

- 

Mass  absorption  coefficient 

for  tissue 

(m2/kfc) 

MUTRHO 

s 

Mass  attentuation  coefficient 

for  air 

(m2/kg) 

MUT 

MS 

Attenuation  coefficient  for  air 

(m~ 1 ) 

DRATE 

8 

Dose  rate  at  one  hour 

(rads/sec) 

TDOSE 

m 

Total  tissue  dose  to  aircrew 

(rem) 

RHOl 

s 

Density  inside  aircraft 

(kg/m3) 

DV 

= 

Dynamic  viscosity 

(Newton- s  ec/m2 ) 

KINV 

8 

Kinematic  viscosity 

(m2/s) 

SOUND 

S 

speed  of  sound 

(m/s) 

DRAGCF 

8 

Slip  correction  factor 

(dimensionless) 

T 

8 

Temperature  at  altitude 

(°K) 

RE 

8 

Reynolds  number 

(dimensionless) 

TERMV 

8 

Terminal  velocity  of  particle 

(m/s) 

RC 

8 

R*Cd 

(dimensionless) 
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APPENDIX  B 


Dose  Program 

This  appendix  contains  the  dose  program  explanation 
and  listing.  The  dose  program  is  designed  to  be  run  either 
interactively  or  by  a  card  deck.  It  requires  five  different 
input  parameters  in  order  to  run.  These  parameters  are: 

(a)  Fission  fraction 

(b)  Yield  in  kilotons 

(c)  Range  oJ  aircraft  to  the  burst  in  meters 

(d)  True  airspeed  velocity  of  the  aircraft  in 
meters  per  second 

(e)  Flight  time  remaining  after  cloud  arrival  in  hours. 

The  program  will  ask  the  user  for  these  parameters  when  they 
are  needed.  The  program  then  calculates  and  displays  the 
aircrafts  cabin  altitude,  the  total  dose  accumulated  in 
rems  for  the  stated  mission,  the  sky-shine  dose  in  rems,  and 
the  activity  density  of  the  outside  air  in  disintegrations 
per  second  per  cubic  meter.  The  information  is  displayed  for 
30  different  altitudes  ranging  from  9150  meters  (approximately 
30,000  feet)  to  305  meters  (approximately  1,000  feet). 

Subroutine  ACT  computes  the  total  outside  activity  at  a 
given  altitude  by  determining  which  particle-size  groups  are 
contributing  to  this  activity.  This  information  is  then 
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passed  on  to  subroutines  SHINE  and  BDOSE  so  that  sky-shine 
and  aircrew  body  dose  can  be  calculated.  Particle  fall 
velocities  and  times  are  calculated  by  subroutine  CLOUD. 
United  States  Standard  Atmospheric  information  is  also 
calculated  by  this  subroutine.  Subroutine  TCHEK  merely 
cross  references  the  time  of  arrival  of  the  aircraft  to  the 
cloud  with  the  distances  that  each  particle-size  group 
distribution  center  has  traveled  in  that  given  time.  The 
subroutine  DENS  obtains  the  correct  air  density  for  the 
corresponding  cabin  altitude. 


PROGRAM  DOSE 

THIS  PROGRAM  CALCULATES  THE  ACTIVITY*  IN  DISINTEGRATIONS  PER 
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APPENDIX  C 


Samp 1 e  Output 


DATA 


INPUT  FISSION  FRACTION 
INPUT  YIELD  (KT) 

INPUT  RANGE  OF  A/C  (M) 

INPUT  A/C  TRUE  AIRSPEED  (M/SI 

INPUT  FLIGHT  TIME  REMAINING  AFTER  TA  IN  HOURS 

#*#*«««#*«##*****#*******«**«****«****##***#»##*#*«* 

TIME  OF  ARRIVAL  <S>  *  3600. 

TIME  OF  ARRIVAL  (HR)  *  1. 

CLOUD  CENTER  HEIGHT  IN  K1LOFEET  »  42.799438 

CLOUD  CENTER  HEIGHT  IN  METERS  =  13045.2687024 

ACPALT  =  903.67449 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
2.359422090052 

SKY-SHINE  DOSE  =  5 . 790639340872 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.000121  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  9150.  IS  4 . 82031 0422S25E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  8.150061420924 

I**************************************************** 

«#******-»******#»«*«*#**•**#**•*******««***»****#»«* 

ACPALT  =  783.14959 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
2.100825493734 

SKY-SHINE  DOSE  =  5.1559756329 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
♦ 00C\255  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  3340.  IS  4 . 455084500 193E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  7.256801126635 

ft*****#*******#*#*********#****#********************. 


****«***«#**#**#***#***«**********#**##*#*******#*** 

ACPALT  =  664.56864 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
1.865349614413 

SKY-SHINE  DOSE  =  4.573056191476 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
,000137  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  8535.  IS  4. 10226901 092E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  6.443405805889 

ft*************************************************** 
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•ft************************************************** 

ACPALT  «  545*  98/69 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  3,  HR  MISSION  IS 
1.534641000955 

SKY-SHINE  DOSE  *  3.766410694189 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.0001426  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  3230.  IS  3, 49S923539154E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 


HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  5.301051695143 
**************************************************** 

**************************************************** 

ACPALT  =  427.40674 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  3.  HR  MISSION  IS 
1.308926458731 

SKY-SHINE  DOSE  =  3.212448129039 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.00014962  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  7925.  IS  3. 092963227236E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  3. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  4. 52137453782 
*****************************************  *********** 


**************************************************** 

ACPALT  *  303.82579 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
1 ♦ 159524984609 

SKY-SHINE  DOSE  =  2. 34577780699 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.000157  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  7620.  IS  2. E38S456 1 S447E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  4.005302791599 

****************************************************  i 

1 

i 

A 
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♦  *m**MI#**t*l*#****«*MH**»**U***»t»*t#*»«t*#»t 

ACPALT  »  190.24484 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
1,042698759571 

SKY-SHINE  DOSE  »  2,559055672582 

MOST  PROMINENT  PARTICLE  SIZE  CROUP  AT  THIS  ALTITUDE  IS 
.0001634  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  7315.  IS  2 . 644205022266E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  3.601754432154 
*«*««***  *****#*#**##****#*#*#*#**#.###**#**  *  #******## 


ft*************************************************** 

ACPALT  =  71.66389 

T$fAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
,.9*41630336135 

iiKY- SHINE  DOSE  =  2.322131893291 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
,00017005  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  7010.  IS  2. 484571 082904E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8, 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  3.268294991904 
**************************************************** 


ft************************* ************************** 

ACPALT  *  0. 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
.8670248018576 

SKY-SHINE  DOSE  =  2. 127905799347 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.00017841  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  6710,  IS  2. 3555739 1873E+9 


TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  2,994930601204 
a*************************************************** 


ft*************************************************** 

ACPALT  =  0. 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  S.  HR  MISSION  IS 
,78784713441 13 

SKY-SHINE  DOSE  =  1,933583079424 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
,00018718  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  6400.  IS  2. 216447541793E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  2,721430213835 
a*************************************************** 

*  *  *  *  ************************************************ 

ACPALT  =  0 . 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
,7301364179016 

SKY-SHINE  DOSE  =  1 . 79 1 94587587 1 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.00019638  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  6100.  IS  2. 124043784637E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  2 . 522032293773 
*****************  * ********************************** 

**************************************************** 

ACPALT  =  0. 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
,6715175027912 

SKY-SHINE  DOSE  =  1 . 648079715241 

MOST  PROMINENT  PARTICLE  SIZE  CROUP  AT  THIS  ALTITUDE  IS 
,00019638  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  5800.  IS  2 . 0 1 95 1 49 1 6699E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  2.319597218032 

ft*************************************************** 


**************************************************** 

ACPALT  =  0 . 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
,622167930959 

SKY-SHINE  DOSE  =  1.526962949178 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
,00020604  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  5490.  IS  1 . 935938497283E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8, 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  2.149130880137 
**************************************************** 
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##***#****#* fr*«***#ft«#***#*#****#ft**«**#*#ft***#*#*#* 

ACPALT  =0, 

TOTAJh  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
.5770797373751 

SKY-SHINE  DOSE  =  1 . 4 1 63047207 1 3 

MOST  PROMINENT  PARTICLE  SIZE  CROUP  AT  THIS  ALTITUDE  IS 
.00021652  DIAMETER 


TOTAL  ACTIVITY  AT  ALTITUDE  5180.  IS  1 . 85736 1 1 43592E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  1.933384458088 
♦a#*#***#*##*****#****#*#******##*******#*********** 


**************************************************** 

ACPALT  =  0. 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
.536663228086 

SKY-SHINE  DOSE  =  1 . 3171 12028276 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.00022862  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  4830.  IS  1 . 7842^6433983E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  1 . 853775256362 
****************************  ***# ******************** 


***************************************************** 

ACPALT  =  0 . 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
,,498442360021  1 

SKY-SHINE  DOSE  =  1.223308014092 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.00022362  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  4570,  IS  1 . 7 1 3239520227E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8, 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  1 . 7217503741 13 

ft*************************************************** 
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ACPALT  =0, 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  S.  HR  MISSION  IS 
*  463522274 1813 

SKY-SHINE  DOSE  =  1 . 137604933437 

MOST  PROMINENT  PARTICLE  SIZE  CROUP  AT  THIS  ALTITUDE  IS 
.00024179  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  4270.  IS  1 . S44942873289E  +  9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  S. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  1 . 601 127257618 
**************************************************** 

**************************************************** 

ACPALT  -  0. 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
. 4305220802764 

SKY-SHINE  DOSE  =  1.056613870104 

MOST  PROMINENT  PARTICLE  SIZE  CROUP  AT  THIS  ALTITUDE  IS 
.00025368  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  3960.  IS  1 . 578724020659E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  1.48713595038 
**************************************************** 


**************************************************** 

ACPALT  =  0 . 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
.4012181864739 

SKY-SHINE  DOSE  =  .9846944447063 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 


.00023368  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  3660.  IS  1 . 518306S87751E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  1.385912631185 
**************************************************** 
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******  ‘  ********************************************* 

ACPALT  =  0. 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
*  3779560595926 

SKY-SHINE  DOSE  =  .9276030966844 

MOST  PROMINENT  PARTICLE  SIZE  CROUP  AT  THIS  ALTITUDE  IS 
.00026829  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  3355.  IS  1 . 4764303SI812E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  1 . 305559 1 56277 
**************************************************** 


**************************************************** 

ACPALT  =  0. 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
.  3525437323577 

SKY-SHINE  DOSE  =  .8652345942124 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.00026829  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  3050.  IS  1 . 421 264792797E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  S. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  1 . 21777832657 
**************************************************** 


**************************************************** 

ACPALT  =0. 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
.3284837776951 

SKY-SHINE  DOSE  =  .8061851680035 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.00028374  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  2745.  IS  1 . 3563630 1 8789E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  1.134668945699 
**************************************************** 


**************************************************** 


ACPALT  =  0. 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8,  HR  MISSION  IS 
.3101831319877 

SKY-SHINE  DOSE  =  . 76 1 2705934 1 56 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.00028374  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  2440,  IS  1 . 330948669312E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  1.071453725403 
**************************************************** 


##***#***************#****#****#*##**###*#*##**#**** 

ACPALT  =  0. 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
,2897135933697 

SKY-SHINE  DOSE  =  , 7 1 1 0329879378 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.0003025  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  2135.  IS  1 . 282050092228E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  1.000746581303 

it*************************************************** 
t*************************************************** 

ACPALT  =  0. 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
.2702162089122 

SKY-SHINE  DOSE  =  .6631813032219 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.00031737  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  1830.  IS  1 . 2329457240SE+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  ,9333975121341 
**************************************************** 

**************************************************** 

ACPALT  =  0 . 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
.2562640813347 

SKY-SHINE  DOSE  =  .6289391303083 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
,00031737  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  1525.  IS  1 . 20537331 7856E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  .8852032116435 
**************************************************** 


♦I#*****#**#***#*#***#*****##***##*##***##***##***#* 

ACPALT  =  0 . 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
.2394179318636 

SKY-SHINE  DOSE  =  , 5875942701853 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.00033835  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  1220.  IS  1 . 160641 163277E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  .8270122020439 
**************************************************** 

♦♦tt**t**«*****t***»**»**«*t*«******««««*»******i*** 

ACPALT  =  0. 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
.2269109068974 

SKY-SHINE  DOSE  =  , 5568987573223 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.00033835  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  915.  IS  1 . 1 33474 1 92062E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  . 7833096642197 

♦  IHHHHHHUt*#*****#*##*******##***##*****#**#**#*#**** 

**************************************************** 

ACPALT  =0. 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
,2116215255932 


SKY-SHINE  DOSE  =  . 5193746137498 

MOST  PROMINENT  PARTICLE  SIZE  CROUP  AT  THIS  ALTITUDE  IS 
.00036071  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  610.  IS  1 . 089030 1 02543E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  .730996139343 
**************************************************** 

a*************************************************** 

ACPALT  =0. 

TOTAL  TISSUE  DOSE  FROM  PARTICLES  INSIDE  A/C  ON  8.  HR  MISSION  IS 
. 1968728070704 

SKY-SHINE  DOSE  =  , 4831773981566 

MOST  PROMINENT  PARTICLE  SIZE  GROUP  AT  THIS  ALTITUDE  IS 
.00036071  DIAMETER 

TOTAL  ACTIVITY  AT  ALTITUDE  305.  IS  1 . 04351 72556 12E+9 
TOTAL  TISSUE  DOSE  ACCUMULATED  IN  REMS  FOR  8. 

HOUR  MISSION  AFTER  FLYING  INTO  CLOUD  IS  .680050205227 
a*************************************************** 


66 


Vita 


Burl  E.  Hickman  was  born  28  August  1945  in  Bryan,  Texas. 
Soon  after,  his  family  relocated  in  Houston,  Texas,  where  he 
completed  high  school  in  1963.  He  entered  Howard  University 
School  of  Engineering  and  Architecture  in  the  mechanical 
engineering  curriculum.  In  his  junior  year,  he  enrolled  in 
the  advanced  ROTC  program  and  received  his  B.S.M.E.  and  USAF 
commission  in  June  1968.  His  first  assignment  was  a  Group 
Mechanical  Engineer  in  a  research  and  development  space 
operations  program  at  Vandenberg  AFB,  California.  During 
this  assignment,  he  earned  a  Master’s  Degree  in  Systems  Man¬ 
agement  from  the  University  of  Southern  California.  In 
January  1973,  he  attended  undergraduate  pilot  training,  and 
upon  graduation  was  assigned  to  Peterson  Field,  Colorado, 
where  he  was  an  Instructor  Pilot  in  the  T-33  aircraft  from 
June  1974  until  May  1976.  He  was  then  transferred  to  K.I. 
Sawyer  AFB,  Michigan.  He  flew  the  F-106  from  May  1976  to  July 
1980,  first  as  a  Squadron  Flight  Commander  and  later  as  the 
Assistant  Operations  Officer.  Major  Hickman  was  assigned  to 
the  Air  Force  Institute  of  Technology's  master's  degree 
program  in  Nuclear  Effects  in  August  1980. 


Permanent  Address:  8330  Lettie  Street 

Houston  TX  77075 


SECURITY  CLASSIFICATION  OF  THIS  PACE  !H»i«i  Data  Enimttd) 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NUMBER 

AFIT/GNE/ PH/82-9 
«.  titlk  (•*>*  Submit) 


2.  OOVT  ACCESSION  NO. 


AIRCREW  ICNIZING  DOSES  FRCN  RADIOACTIVE 
DUST  CLOUD  G0JERATH)  BY  NUCLEAR  BURST 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1.  RECIPIENT’S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  A  PERIOD  COVERED 


MS  Thesis 


S.  PERFORMING  ORG.  REPORT  NUMBER 


AUTHOR!  R) 


I.  CONTRACT  OR  GRANT  NUMBER!*) 


Burl  E.  Hickman 
Major  USAF 

S.  PERFORMING  organization  name  and  ADDRESS  10.  PROGRAM  element,  project,  task 

AREA  A  WORK  UNIT  NUMBERS 

Air  Force  Institute  of  Technology  (AFIT-Bl) 

Wright-Patterson  Air  Force  Base,  Ohio  45433 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  U.  REPORT  DATE 

March  1982 _ 

IS.  NUMBER  OF  PAGES 

_ _ 1A _ 

14.  MONITORING  AGENCY  NAME  &  ADDRESS!// di/Zarant  /ram  Controlflnf  Ottlcm)  IS.  SECURITY  CLASS,  (ot  thtt  rtport) 


ISa.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 

is.  DISTRIBUTION  STATEMENT  (o /  Ihlt  Report) 


Approved  for  Public  Release;  Distribution  Unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  ebatrmct  entered  in  Block  20,  it  different  from  Report) 


1  b  APR  i««2 


I  IB.  SUPPLEMENTARY  NOTES 


APPROVED  r 


«rR  190-1/. 


UCufi  i or  i 0i i  ■-..ij 

Professional  Development 
Air  Force  Institute  of  Technology  (ATC) 
Wright-Patterson  AFB,  OH  45433 


19.  KEY  WORDS  fConr/nu*  on  reverse  side  it  neceeeery  end  identity  by  block  number ) 

Nuclear  Radiation 

Dust 

Crew  Survivability /Vulnerability 

Radioactive  Dust  Clouds 

Cloud  Penetration _ 

20.  ABSTRACT  (Confirm#  on  reverae  aide  it  neceeeery  end  Identity  by  block  numbmt) 

This  report  will  evaluate  the  threat  of  radioactive  fallout  to  which  aircrew 
members  will  be  ejqjosed  when  flying  through  a  descending  fallout  cloud.  A 
computer  program  is  developed  for  calculating  the  ionizing  dose  rate  of  a 
radioactive  dust  cloud  as  a  function  of  time,  and  also  the  dose  that  an  aircrew 
receives  when  flying  through  the  respective  cloud.  A  cloud  model  that  is 
patterned  after  the  AFIT  fallout  smearing  code  was  developed.  A  comparison  is 
made  between  the  activities  at  various  altitudes  from  305  meters  tO  9150  meters 
_ f Continued  on  Reversed _ 

DD  1  JAN *71  1473  EDITION  OF  I  NOV  45  IS  OBSOLETE 


■«'*■**-  vwjv teLmtsBkr. 


SECURITY  CLASSIFICATION  OF  THIS  RAGE  (When  Dete  Bntoted) 


_ UNCLASSIFIED _ 

WCUWITV  CLASSIFICATION  OF  THU  PAOtflWiaw  Data  InUwW) _ _ 

BLOCK  20:  ABSTRACT  (Continued) 

to  provide  information  for  possible  re-direction  of  flight.  The  external 
ionizing  dose  to  the  aircrew  is  computed  by  the  new  code  considering  the 
cloud  size,  the  aircraft's  transit  time  through  the  cloud,  and  the  ingestion 
rate  of  radioactive  particles  into  the  aircraft’s  cabin.  Information  is 
also  provided  to  indicate  the  method  by  which  doses  can  be  computed  from  a 
cloud  of  multiple  bursts.  The  results  demonstrate  that  total  dose  to  each 
aircrew  member  is  approximately  8  rems  aft«r  flying  through  a  fallout  cloud 
one  hour  after  cloud  stabilization  of  a  1  hi  burst,  with  the  mission  con¬ 
tinuing  for  eight  hours  subsequent  to  the  cloud  transit. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfWian  Data  EntararfJ 


